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We consider the process of Higgs boson pair production at a future proton collider with centre- 
of-mass energy of 100 TeV, focusing on rare final states that include a bottom-anti-bottom quark 
pair and multiple isolated leptons: hh —>■ {bb) +n£ X, n = {2,4}, X = —}. We construct 

experimental search strategies for observing the process through these channels and make suggestions 
on the desired requirements for the detector design of the future collider. 


I. INTRODUCTION 

In the ‘aftermath’ of the Large Hadron Collider’s 
(LHC) ‘Run 1’ at 7 and 8 TeV, we are left with the 
grim possibility that striking new effects will be beyond 
the reach of the collider, perhaps even in ‘Run 2’ at 
13 TeV. The LHC will continue to collect data into a 
high-luminosity ‘Run 3’ phase, possibly at 14 TeV, which 
will allow us to perform precision tests on the properties 
of the Higgs boson and the underlying mechanism of elec- 
troweak symmetry breaking (EWSB). Taking a rather 
pessimistic point of view when it comes to beyond-the- 
Standard Model physics, a future circular hadron collider 
(FCC-hh), at a potential proton-proton centre-of-mass 
energy of 100 TeV, will resume this task, zooming into 
the properties of the Higgs boson. An important pro¬ 
cess that has received considerable attention in the re¬ 
cent years, as part of this effort, is that of Higgs boson 
pair production. Examining this process, several aspects 
of EWSB can be probed: e.g. the consistency of the self¬ 
coupling of the Higgs boson with the Standard Model 
expectation [IHH], whether it comes as part of a doublet 
or is affected by new physics effects as well as 

effects of higher-dimensional operators [451IT7] . 

The hh process is currently under investigation at the 
LHC, both by phenomenologists and experimentalists. 
The ATLAS and CMS collaborations have already pre¬ 
sented, as a “warm-up” exercise, results for resonant and 
non-resonant production of hh using ‘Run 1’ results |48l - 
Eo]. Detection of this process remains, however, a se¬ 
rious challenge, even at the end of the high-luminosity 
run of the LHC, where 3 ab“^ of integrated luminosity 
will be collected. The question examined in this article is 
whether a 100 TeV proton-proton collider can yield sig¬ 
nificant contributions to the investigation of this process. 

At present, it is clear that the task will not be trivial, 
even at the FCC-hh, but initial results encourage op¬ 
timism. Several studies examining the hh —>■ (bb){'yj) 
channel, with diverse assumptions on detector perfor¬ 
mance, indicate that it would provide a clear signal at 
the end of a 3 ab“^ run of the FCC-hh [321 [ST]. This 


* Electronic address: apapaefs@cern.ch 


process falls into the ‘rare but clean’ category, allowing 
good rejection of relatively manageable backgrounds with 
high signal efficiency, on account of the prospect of excel¬ 
lent resolution in the photon momentum measurement. 
An obvious question is whether the increase of energy 
opens up access to other such channels, previously un¬ 
available at the LHC due to the small total cross sec¬ 
tion. In the present article we consider a simple set of 
such channels, where one Higgs boson is allowed to de¬ 
cay via h ^ bb and the other to either a pair of gauge 
bosons h —>■ ZZ,Zj,W~^W~, with subsequent decay to 
leptons, or directly to leptons h —>• (T+r“)/(/r+^“). We 
will not consider hadronic decays of the r lepton, as these 
will require reconsideration of tagging algorithms versus 
large QCD backgrounds that will rely crucially on the, 
presently unknown, detector performance of the future 
collider.* 

The final states considered, their branching ratios and 
cross sections at 14 TeV are shown in Table HI We also 
give the relevant numbers for the hh —>■ (&5)(77) process, 
even though it has not been considered here. Evidently, 
of the multilepton channels examined here, all but the 
hh —)■ {hb){W'^W~) channel provide a negligible num¬ 
ber of events at the LHC at 14 TeV. These calculations, 
and the rest of the signal cross sections in this article, are 
based on the EFT-approximate total NNLO cross section 
for Higgs boson pair production of (t( 14 TeV) = 40.2 fb 
and cr(100 TeV) = 1638 fb [53]. At the FCC-hh, the 
factor ^ 40 increase in cross section, with respect to the 
LHC at 14 TeV, allows for exploration of these rare chan¬ 
nels, given that they provide a non-negligible inclusive 
event sample. In what follows and the rest of this article, 
we define £ = (e, pt}. 

It is important to mention that the theoretical uncer¬ 
tainties on both the signal cross section and differential 
distributions are at present rather large. These include 
0(10%) due to scale uncertainty, 0(10%) due to parton 
density function uncertainty [23 ED] and an estimated 


*We note also the recent study at the LHC at 14 and the FCC- 
hh at 100 TeV of the three-lepton final state arising from hh 
{W^W-){W^W-) that appeared in m- 
^For further details on the cross section calculations see also |421I54F 

Eg. 
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0 ( 10 %) due to the unknown finite top mass effects at 
next-to-leading order in QCD [42]. At present, the ex¬ 
act electroweak corrections at next-to-leading order also 
remain uncertain, estimated to be 0(5%) for the case of 
single Higgs boson production at the LHC m- We ex¬ 
pect that higher-order predictions for the signal as well as 
the backgrounds, both QCD and electroweak, at next-to- 
leading order and beyond, will be substantially improved 
in time for the construction and data taking of a future 
hadron collider, which is foreseen to take place in 20 - 
30 years, and we will neglect them in this initial study 
for the sake of simplicity. Moreover, the current esti¬ 
mate of 0(30%) on the total theoretical uncertainty are 
not expected to substantially alter the conclusions of the 
analyses. 

The paper is organised as follows: we first give details 
of the simulation of detector effects and the generation of 
Monte Carlo event samples used for the analyses. Subse¬ 
quently, we go through each final state, highlighting the 
important attributes and potential for measurement at 
the FCC-hh. The main results, in Tables 


and VH appear as expected number of events at 3 ab 
of integrated luminosity and can be trivially translated 
into cross sections or events at higher luminosities. We 
conclude by summarizing and making general recommen¬ 
dations on desired features of the FCC-hh detectors. 


II. EVENT GENERATION AND DETECTOR 
SIMULATION 


A. Detector simulation 

In what follows, we consider all particles within a pseu¬ 
dorapidity of \r]\ < 5 and pr > 500 MeV. Following |5T] . 
we smear the momenta of all reconstructed objects ac¬ 
cording to [62] for jets and muons and for electrons ac¬ 
cording to [53] . We also apply the relevant reconstruction 
efficiencies according to [52] for jets and muons and [53] 
for electrons. We simulate 6 -jet tagging by looking for 
jets containing H-hadrons and applying a flat b-tagging 
efficiency of 70% and a mis-tag rate of 1% for light-flavour 
jets. We do not consider c-jets. We do not apply any 
smearing to the missing transverse energy. 

We reconstruct jets using the anti-Zcj algorithm avail¬ 
able in the Fast Jet package [HU [55], with a radius 
parameter of i? = 0.4. We only consider jets with 
Pt > 40 GeV within \r]\ < 3 in our analysis. The 
jet-to-lepton mis-identification probability is taken to be 
= 0.0048 X e-o o35PT,/GeV^ in m- We apply a 
transverse momentum cut on all leptons of px > 20 GeV 
and require them to lie within a pseudorapidity \r]\ < 2.5. 
We also consider the mis-tagging of two bottom quarks 
with a flat probability of 1 % for each mis-tag, correspond¬ 
ing to a 6 -jet identification rate of 70% and demand that 
they lie within I 77 I < 2.5. We demand all leptons to be 
isolated, an isolated lepton having ^ ■ px^i less than 15% 
of its transverse momentum in a cone of Ai? = 0.2 around 


it. 

Since the design of the detectors for FGG-hh is still 
under development, we consider in conjunction to the 
above, which we call ‘LHC parametrization’ of the detec¬ 
tor effects, an alternative ‘ideal’ parametrization. This 
is obtained by setting all efficiencies to 100%, within 
the given considered acceptance regions for jets and lep¬ 
tons, and by removing all momentum smearing effects. 
Table Hi summarizes the differences between the two 
parametrizations. The mis-tagging rates for 6-jets, lep¬ 
tons and photon were kept identical in both parametriza¬ 
tions. We note, however, that these were not found to 
be particularly important for the channels considered 
here, provided that they remain at the levels that can 
be achieved at the high-luminosity LHC. 


B. Event generation 

We generate the signal at leading order using the 
Herwig++ event generator [661169] interfaced with the 
OpenLoops package for the one-loop amplitudes [iniiTg.* 
The backgrounds have been generated with the MadGraph 
5/aMC0NLO package [7TH73] . at next-to-leading order 
(NLO) in QCD, except for the case of the ti background, 
which was generated at leading order and merged with 
the Herwig++ parton shower using the MLM algorithm, 
including parton matrix elements. For the latter, we 
normalized the cross section to the total NLO cross sec¬ 
tion. All simulations include modelling of hadronization 
as well as the underlying event through multiple parton 
interactions, as they are made available in Herwig++. No 
simulation of additional interacting protons (pile-up) is 
included in this study. The CTlOnlo parton density func¬ 
tion m set was used for all parts of the simulation. 

To facilitate the Monte Carlo event generation, we re¬ 
duce the rather large tt cross section at 100 TeV by apply¬ 
ing the following generation-level cuts to the final state 
objects (£+6r'^)(^'“6i/^/):§ 

PT,b > 40 GeV, pT,i > 30 GeV, 

\m\ < 2.5, 

0.1 < A(6,6) < 2.0,0.1 < A(£, Q < 2.0. (1) 

We also apply generation-level cuts to 6-quarks in other 
backgrounds. These are indicated on the corresponding 
tables. 


^Note that 0(10%) systematic uncertainties can be introduced by 
using the showered LO sample instead of a sample that includes 
higher-order real emission matrix elements. An extensive discus¬ 
sion of the parton shower uncertainties in Higgs boson pair produc¬ 
tion can be found in m- For the purposes of this initial study, it 
is sufficient to use LO hh production. 

^We have checked that a full MC@NLO tt sample produces similar 
efficiencies before running out of statistics. 
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channel 

BR 

cr(14 TeV) (fb) 

a(100 TeV) (fb) 

hh -4 {bb){ZZ) {bb){e+£-£+£-) 

0.016% 

0.006 

0.26 

hh ->■ {bb){Z-f) {bb){£+£--f) 

0.013% 

0.005 

0.21 

hh lbb)(W+W-) {bb){£+£'- + $) 

1.658% 

0.667 

27.16 

hh -4 (55)(r+r") ^ {bb){£+£'- + $) 

0.893% 

0.360 

14.63 

hh —>■ {bb){p^ p,~) 

0.025% 

0.010 

0.42 

hh —>• (55) ( 77 ) 

0.263% 

0.106 

4.31 


TABLE I: Higgs boson pair production rare final state branching ratios (BR) and cross sections {a, in femtobarn) that are 
considered in the present article. Decays to r leptons that subsequently decay to electrons or muons are included. The 
hh —>■ {bb) {')')) channel is not considered here and given for comparison. 


property 

LHC param. 

ideal param. 

Pb —tag 

70% 

80% 

e(i) 

Ref. |62J 

100 % 

e(e) 

Ref. [63] 

100 % 

e(h) 

Ref [62] 

100 % 

O' O') 

Ref [62] 

0 

o-(e) 

Ref [63| 

0 


Ref. 

0 


TABLE II: The differences between the two detector 
parametrizations we consider. One is effectively an ‘ideal’ de¬ 
tector, with perfect efficiency within the considered pseudora¬ 
pidity ranges for each type of object, whereas the other is an 
LHC-like parametrisation, with equivalent assumptions as for 
the high-luminosity LHC. The 6 -tagging probability, Pb-tag, 
the efficiency for jets, electrons and muons (e(i/e/p)) and the 
corresponding standard deviation of the smearing Gaussian 
are given {a{j/e/fi)). 


III. ANALYSIS 

A. hh^{bb){4£) 

Since the h —>■ ZZ —>■ M channel has played an im¬ 
portant role in the discovery of the Higgs boson, it is 
reasonable to ask whether the four-lepton final state, in 
association with h (56), could be employed equiva¬ 
lently in the di-Higgs channel. Unfortunately, at the 
LHC at 14 TeV, the SM cross section in this channel 
is hopelessly low: cr((55) -|- 4£)i4 xeV ~ 6 x 10“^ fb, giv¬ 
ing an expected number of ~ 20 events at the end of 
the high-luminosity run at 3 ab“^. Hence, even consid¬ 
ering just triggering and basic acceptance cuts, one can 
conclude that this channel will never be observed at the 
LHC. At a 100 TeV collider, the cross section increases 
to about 0.26 fb, leading to ^ 780 events at 3 ab“^. 
Evidently, the channel is still challenging, even at this 
future collider. Nevertheless, the final state is easier to 
reconstruct than others, and one should consider whether 
significance could in principle be obtained. Particularly 
interesting is the scenario of an integrated luminosity of 
30 ab“^, where one would have an initial sample of 7800 
events. 

The backgrounds relevant to this process are listed in 


Table m Here, we only consider mis-tagging of a sin¬ 
gle lepton, with the the dominant process in this case 
being W^Zh. Evidently, there are 6 objects relevant 
to the hard process: two 5-jets and 4 leptons. We de¬ 
mand pairs of leptons of opposite charge and same flavour 
as well as two identified 5-jets. As a simulation of a 
possible 4-lepton trigger, we ask for the following stag¬ 
gered lepton cuts: PT,f{i_ 2 , 3 , 4 } > {35,30,25,20} GeV. 
Since the signal in this case is not expected to possess 
a large amount of missing transverse energy, we impose 
< 100 GeV. It was observed that the distance be¬ 
tween all leptons in the hh signal is substantially smaller 
than in most of the background processes, and hence 
a cut of Ai?(.^i,£j) < 1.0, with j = {2,3,4}, is im¬ 
posed. We apply cuts tailored to rejecting events with 
two on-shell Z bosons: if there are two combinations of 
same-flavour opposite-sign leptons that have an invariant 
mass in mi+i- € (80,100) GeV, we reject the event. We 
also demand that no single pair of same-flavour opposite- 
sign leptons possesses a mass above 120 GeV. The final 
reconstructed observables are the invariant mass of the 
four-leptons, the invariant mass of the 5-jet pair and the 
invariant mass of all six reconstructed objects. These ob¬ 
servables are shown for the hh signal and the two signifi¬ 
cant backgrounds, tth and tiZ, after the aforementioned 
cuts in Fig. To obtain the final result we further im¬ 
pose the following cuts: 

Mbb e (100,150) GeV, e (110,140) GeV, (2) 

and no cut on the invariant mass of all the reconstructed 
objects. 

It is important to emphasize at this point a crucial 
element of this analysis, namely the minimal cuts that 
should be applied each of the four leptons in the final 
state to avoid excessive signal rejection. In particular, 
for the fourth, softest, lepton in the final state, about 
65% of the events fall in 20 GeV < pt/a < 30 GeV. Not 
including this bin in pT,i 4 would automatically reduce the 
signal by more than a factor of two. This is demonstrated 
in Fig. where the ordered transverse momenta of the 
four leptons are shown. The ordered transverse momenta 
of the 5-jets are shown for completeness. 
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channel 

cr(100 TeV) (fb) 

Ns ab-i(cuts, ideal) 

^3ab-i(cuts, LHC) 

hh-> {bb){£+£-£+£-) 

0.26 

13.01°:^ 


tth^ {£+bi^e)i£'~bPp){2£) 

ttZ {£+but)\£'-bvp){2£) 

Zh^ (66)(4.f) 

ZZZ ^ (4£)(66) 

bbh — >• bb{4£), pT,b > 15 GeV 

ZZh^(4f)(66) 

193.6 

256.7 

2.29 

0.53 

0.26 

0.12 

30.4+il 

6.6+?® 

0(10-1) 

0(10-1) 

0(1) 

0(10-®) 

lO.OlK;? 

^•^-0.7 

0(10-1) 

0(10-1) 

0(10-1) 

0(10-®) 

ZZ — >■ (dfj-l-mis-tagged 66 

781.4 

0(10-^) 

0(10-®) 

hZ — >■ (4^)-|-mis-tagged 66 

68.2 

0(10-®) 

0(10-®) 

W ± ZZ — > {£vi){£'^ £~)(bb) -\- mis-tagged £ 

7.5 

0(10-®) 

0(10-®) 

W^Zh — >■ {£i'e){£^£~){bb) -I- mis-tagged £ 

1.4 

0(10-®) 

0(10-®) 


TABLE III: The signal and backgrounds that have been considered for the {bb){£^£~£ ^£ ~) channel are shown in the second 
column. All the background cross sections have been calculated at NLO using the MadGraph5/aMC@NL0 package. We show the 
resulting number of expected events, ab-ii s-f s-n integrated luminosity of 3 ab“^ for the two detector parametrizations in 
the third and fourth columns. We show In-equivalent errors, derived according to the Poisson distribution. 



FIG. 1: The final reconstructed observables for the (66)(4£) final state, for the LHC parametrization: the invariant mass of 
the four-leptons, the invariant mass of the 6-jet pair and the invariant mass of all six reconstructed objects. In the former two, 
the signal has been multiplied by a factor of 10 for clarity. The invariant mass of all objects is given after cuts on the prior 
two observables have been imposed. The Monte Carlo statistical errors are shown for each data sample to help distinguish 
statistical fluctuations from physical features. 




(a) (b) 

FIG. 2: The transverse momenta of the four leptons in the (66) (4£) final state prior to any cuts are shown on the left, for 
the hh signal. The leptons are ordered in pr with pT,ei > PT,ij for j > i. The right-hand figure shows the ordered transverse 
momenta of the 6-jets for completeness. 
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B. hh^ {bb){i+r-i) 

The hh (bh){t^l~^) channel contains two leptons 
from the on-shell Z decay and a hard photon, and pos¬ 
sesses a cross section only slightly lower than hh 
{bb){4:£). The backgrounds are, however, substantially 
larger in this case. Here we only include the most signif¬ 
icant irreducible ones, coming from bbZj, Uj, and hZ"f, 
as well as the dominant reducible ones, where a photon 
is mis-tagged in bbZ or tt production. 

In the analysis of this channel we require two leptons 
of the same flavour with 2 } > {40,35} GeV, two 
anti-kt R = 0.4 6-jets with PT,b{i 2 } > {60,40} GeV, 
< 80 GeV and a photon with pT,'y > 40 GeV. 
No isolation requirements are imposed on the photon. 
We ask for Ai?(t!i,t' 2 ) < 1.8, Ai?(£i, 7 ) < 1.5 and 
0.5 < AR{bi,b 2 ) < 2.0. We construct the invari¬ 
ant mass of the 6-jet pair and the invariant mass of 
the two-lepton and photon system and impose the cuts: 
100 < < 150 GeV and 100 < < 150 GeV. 

Even after these cuts, the bbZj background dominates 
the resulting sample, resulting to a signal-to-background 
ratio of 0{2 — 3%) with only 0(10) signal events at 
3 ab“^ of integrated luminosity. Therefore, this chan¬ 
nel is not expected to provide significant information at 
the 100 TeV pp collider, unless a signihcant alteration of 
the hh channel is manifest due to new physics effects. 


c. hh^ {bb){i+r){+$T) 

Another interesting channel to consider is the final 
state that includes a 66 and two oppositely-charged lep¬ 
tons. This channel receives signal contributions from 
three different hh decay modes. The largest contribution 
comes from hh {hb){W'^W~) with the Ws decaying 
(either directly, or indirectly through taus) to electrons 
or muons. The second-largest contribution to this chan¬ 
nel comes from hh —>• (66)(t+t“), with both taus decay¬ 
ing to electrons or muons. Both of these channels will 
include final-state neutrinos, and hence will be associ¬ 
ated with large missing energy. The specific final state 
has been considered at the LHG at 14 TeV, e.g. com¬ 
ing from {bb)(W+W-) in. [ 4 ] or included implicitly as 
part of the channel (66)(t+t“) in [5H5lll0j. The smallest 
contribution comes from hh —>■ (bb){pApr), i.e. through 
direct decay of one of the Higgs bosons to muons. This 
has been considered in [2] at a 200 TeV proton collider, 
which was envisioned to have a integrated luminosity of 
either 600 fb“^ or 1200 fb“^. The channel was shown to 
be able to provide some information on the hh process 
at a higher-energy collider, and hence we re-examine it 
here at a 100 TeV collider, allowing for the possibility of 
collecting a higher integrated luminosity sample. 

Due to the different origin of the leptons in the three 
different signal processes, the kinematical details vary 
substantially. As already mentioned, in {bb){W~^W~) 


and (66)(r“''T“), we expect large missing energy. The t 
leptons in (66)(T+r“) are light compared to the Higgs 
boson, and hence the leptons and neutrinos in their de¬ 
cays are expected to be collimated. On the other hand, in 
{bb){W^W~), both Ws are heavy, one being most of the 
time on-shell {Mw ~ 80.4 GeV) and the other off-shell 
with Mw* peaking at ~ 40 GeV. This implies that the 
distribution of the leptons and neutrinos will differ con¬ 
siderably from the (66)(r+T“) case. This will be reflected 
in the analysis efficiency for the two processes: the fact 
that the leptons in boosted r decays are collimated with 
the associated neutrinos makes the h r+T“ decay eas¬ 
ier to separate from the backgrounds. The (66)(^+^“) 
final state is expected to have little missing energy, possi¬ 
bly only due to H-hadron decays, and the two muons are 
expected to reconstruct the Higgs boson. To account 
for all of these properties, we construct two separate 
signal regions. One aims to capture events containing 
rather large missing energy, targeting the (66)(W'‘W“) 
and (66)(r“'"r“) channels, whereas the second is aimed 
towards events with minimal missing energy that are ex¬ 
pected to characterize the {bb){pApT) channel. The same 
set of observables is considered, with substantial varia¬ 
tions of cuts. We consider events with two isolated lep¬ 
tons, with isolation criteria as in the previous sections. 
We construct their transverse momentum, 
distance between them, AR{(.l,^ 2 ) and their invariant 
mass Mu. We ask for two 6-jets, for which we construct 
the transverse momenta, distance and invariant mass: 
PT,f){i 2 }) AR{bi,b 2 ), Mbb- We also consider the missing 
transverse energy, and the invariant mass of all the 
reconstructed objects, Mbbu- Moreover, we consider a 
further observable, M^eco.j constructed by assuming that 
the missing energy arising from neutrinos in the decays 
of the T leptons is collinear to the observed leptons: 

Mreco. = [Pbi + Pb2 + (1 + fl)Pli + (1 + /2)W2]^ i (3) 

where pb^ , Pi^ are the observed momenta of the f-th 6-jet 
and i-th lepton and /i _2 are constants of proportional¬ 
ity between the neutrino and lepton momenta from the 
decay of the two r leptons: Pi^i — fiPh- These can be 
calculated from the observed missing transverse energy 
by inverting the missing transverse momentum balance 
relation Lf = where L is the matrix = Pe ^ for 
which the superscript denotes the component of the i-th 
lepton momentum, j = {x, y} and E and f are the vec¬ 
tors $ = and f = (/i,/ 2 )- We consider this 

observable for the {bb){W^W~) signal sample as well, 
even though the collinearity approximation fails, as it is 
still expected to be correlated with the invariant mass 
of the Higgs boson pair. One may also dehne the re¬ 
constructed Higgs boson mass from the reconstructed r 
lepton momenta: Mh.reco, = [(1 + .fi)pii + (1 + f 2 )pi 2 ] ■ 

We call the signal regions SR^ and SR^, corresponding 
to the hh decay modes {bb){W^W~), (66)(r+T“) and 
{bh){pApL~) respectively. Table [v| shows the cuts chosen 
for these observables for these three signal regions. Note 
that the Mt 2 observable can also be constructed for the 


6 


channel 

(t( 100 TeV) (fb) 

N 3 ab-i(cuts, ideal) 

A 3 ,b-i(cuts, LHC) 

hh^ {bb)(e+r-/) 

0.21 

ii.of:);;® 

7 

' •'^-0.2 


bbZ7 bb{£+£-)'r,pT,b > 30 GeV 26.00x10^ 266l™“ 203^^ 

ttj ^ {L+bi^LOiL-buL)^ 7.94x10® TSias 

hZ7 ^ {bh){i+r)-i _L72_3t?2l® 


bbZ —>• bb{t^l )-|-mis-tagged 7 , pT,b > 30 GeV 

107.36x10® 

20 lf 

2 U- 

tt —> (^^6i'^)(/“67'^/)-|-mis-tagged 7 , cuts as in Eq.[^ 

25.08x10® 

14+1 

lOlj 


TABLE IV: The signal and backgrounds that have been considered in the {bh){t^l~^) channel. All the background cross 
sections are given at NLO using the MadGraph5/aMC@NLD package. The tt process with the mis-tagged 7 was generated at 
tree-level, merged to one jet via the MLM method, and normalized to the NLO cross section. We show Icr-equivalent errors, 
derived according to the Poisson distribution. 


rejection of the tt background, but we do not take this 
approach here [ 3 . 

The backgrounds considered for the 
final state include the irreducible ones coming from ti, 
with subsequent semi-leptonic decays of both top quarks, 
those from bbZ with decays of the Z boson to leptons, 
those from bbh with subsequent decays of the Higgs bo¬ 
son to two leptons, as well as the resonant hZ and ZZ 
backgrounds. We also consider here the mis-tagging of a 
jet to a single lepton through the bbW^ channel and the 
mis-tagging of bb in the -l-jets background, which 
was considered using parton at NLO. As before, 

we do not consider mis-identification of c-jets to 6-jets. 


We show the resulting events after analysis for the sig¬ 
nal region SR^ in Table VI Both the “LHC” and the 
“ideal” parametrizations are shown. Evidently, since we 
are using the same set of cuts for both scenaria, the 
“ideal” parametrization does not necessarily provide a 
substantial improvement to the signal efficiency. This ef¬ 
fect is observed in particular for the {bh){W^W~) sam¬ 
ple. Moreover, it is evident that this channel could be ad¬ 
dressed by employing more advanced statistical methods: 
cuts could be devised separately for the two sub-channels 
{bb){W~^W~) and (66)(T+r“) and then combined taking 
into account the correlations and cross-contamination.^ 
For example, this could be done via appropriate cuts on 
Afh.reco.! which was found, as expected, to peak at the 
Higgs boson mass for the {bb){T'^T~) sample and around 
50 GeV for the {bb){W~^W~) sample. This study is be¬ 
yond the scope of this initial investigation. 


Due to the fact that both signal and backgrounds pos¬ 
sess larger cross sections, harder cuts are imposed in this 
channel than in the ( 66 ) (4€) and ( 66 ) (Z 7 ) channels. This 
involves cutting on variables that are expected to have a 
high degree of correlation with the invariant mass of the 
Higgs boson pair, which implies that this channel would 
be less sensitive to variations of the self-coupling than, 
for example, the hh —>■ ( 66 ) (4t'), where no such cuts are 


^Another improvement would involve separating the (fefjjfr+r”) 
final states into same-flavour and different-flavour. The author 
thanks Fabio Maltoni for useful discussion on this. 


imposed. Nevertheless, insofar as Standard Model-like 
hh production is concerned, this process is expected to 
provide important information, contributing to the de¬ 
tection of this channel at a 100 TeV collider, with signal- 
to-background ratio of ~ 0.1 and large statistical signifi¬ 
cance at 3 ab“^. 

On the other hand, the situation for the ( 66 )(/i+/i“) 
after the SR^ cuts are applied is rather bleak: at 3 ab“^ 
of integrated luminosity, only a handful of events are ex¬ 
pected with the “LHC” detector parametrization with a 
few hundred background events, even after hard trans¬ 
verse momentum cuts on the muons and a tight mass 
window on the di-muon invariant mass around the Higgs 
boson mass. Because of the latter cut, turning to the 
“ideal” situation improves the signal efficiency substan¬ 
tially, since the smearing of the muon momenta is absent. 
Despite this, only 0(10) events would be obtained at 3 
ab“® with a similar number of background events as for 
the ‘LHC’ parametrization. Hence, barring any signifi¬ 
cant enhancements of the rate due to new physics, the 
{bb){fi^ fi~) contribution to the hh (bb){(.'^i~) final 
state is not expected to provide significant information. 


IV. DISCUSSION AND CONCLUSIONS 

In this letter we have considered three rare, and poten¬ 
tially clean, final states coming from standard model-like 
Higgs boson pair production at a future 100 TeV proton- 
proton collider. These processes are made available for 
investigation due to the fact that the total Higgs boson 
pair production cross section increases by a factor of ~ 40 
over that expected at the LHC. 

For the hh —>■ {bh){ZZ) channel, where both Z bosons 
decay to leptons, resulting in ( 66 )(4£), it was shown that 
a few events could be obtained at 3 ab“®, versus 0(10) 
background events. A crucial aspect of this analysis, 
that should be taken into account in detector design, 
is being able to observe four-lepton final states with at 
least one lepton possessing transverse momentum down 
to - 20 GeV. 

The hh —>■ ( 66 ) (Z 7 ) channel, with the Z boson decay¬ 
ing to leptons, was also considered briefly and was found 
to be of negligible importance even at 30 ab“®. This is 
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observable 

SRg 


SR^ 

f^T 

> 100 GeV 


< 40 GeV 

PT,ti 

> 60 GeV 


> 90 GeV 

PT,12 

> 55 GeV 


> 60 GeV 

ARi£i,e2) 

< 0.9 


G ( 1 . 0 , 1 . 8 ) 


€ (50,80) GeV 


G (120,130) GeV 

Ptm 

> 90 GeV 


> 90 GeV 

Ptm 

> 80 GeV 


> 80 GeV 

Ai?( 6 i,b 2 ) 

G (0.5,1.3) 


G (0.5,1.5) 

Mbbu 

> 350 GeV 


> 350 GeV 

Mbb 

e (110,140) GeV 


G (110,140) GeV 

^^reco. 

> 600 GeV 


none 

TABLE V: The cuts for the three signal regions constructed in the analysis of the hh —>■ (bb)(£^£ 

) channel. The signal 

regions SR^ and SR^, correspond to the hh decay modes {bb){W'^W ), {bb){' 

r'^T ) and {bb){pL^p ) respectively. 

channel 

a(100 TeV) (fb) 

^3 ab-i(cuts, ideal) 

iV 3 ,b-i(cuts, LHG) 

hh {bb){W+W-) {bb){£'+vtR-Pt) 

27.16 

20.9ti'^ 

19.97i:^ 

hh —> (fe 6 )(r"'"r“) —>■ {bb){£''^viiVr£~P iVt) 

14.63 

38.511^ 

24.31®;® 

tt —> {£'^bvi){£'~bvii), cuts as in Eq. 

25.08x10^ 

34.3tf4^ 

15.8lf8® 

bbZ ^ bb{£+£-), pT,b > 30 GeV 

107.36x10^ 

257.9l74®6^ 

493.7l“t.4 

ZZ ^ bb{£+£-) 

356.0 

0 ( 10 -®) 

0 ( 10 -®) 

hZ bb{£+£-) 

99.79 

49.8+1^ 

40.411® 

bbh bb{£+£-), pT,b > 30 GeV 

26.81 

0 ( 1 ) 

0 ( 1 ) 

fabW^*^ —>■ bb{£^i'i), pT,b > 30 GeV -|-mis-tagged £ 

1032.6 

0 ( 10 -^) 

0 ( 10 -®) 

£'^t~ -1-jets —^ {V'£~) - 1 - mis-tagged bb 

2.14x10® 

0 ( 10 -®) 

0 ( 10 -®) 

TABLE VI: The signal and backgrounds that have been considered in the {bb){£'^£~ + $) channel coming from h —>■ 

All the background cross sections were calculated at NLO using the MadGraphS/aMCSNLO package. The tt process was generated 

at tree-level, merged to one jet via the MLM method, and normalized to the NLO cross section. We show la-equivalent errors. 

derived according to the Poisson distribution. 




channel 

a(100 TeV) (fb) 

^3 ab-i(cuts, ideal) 

7^3 ab-i(cuts, LHC) 

hh —> 

0.42 

8 . 6 :i:i 

1 Q+U.b 

-*-•^-0.3 

tt —^ {£'^bvi){£'~bvii), cuts as in Eq. 

25.08x10® 

48.07;i®'o" 

15.8lfg^ 

bbZ ^ bb{£+£-), pT,b > 30 GeV 

107.36x10® 

< 73.5 

49.4l®4®®4" 

ZZ bb{£+£-) 

356.0 

0 ( 10 -®) 

0 ( 10 -®) 

hZ ^ bb{£+£-) 

99.79 

0 ( 10 -®) 

2.5l°;I 

bbh bb{£+£-), pT,b > 30 GeV 

26.81 

0 ( 1 ) 

0 ( 1 ) 

bbW‘’= —>• bb{£^Vf_), pT,b > 30 GeV -|-mis-tagged £ 

1032.6 

0 ( 10 -®) 

0 ( 10 -®) 

l~ -1-jets —^ (£'^£~^ -1- mis-tagged bb 

2.14x10® 

0 ( 10 -®) 

0 ( 10 -®) 


TABLE VII: The signal and backgrounds that have been considered in the jj,~) channel. All the background cross 

sections were calculated at NLO using the MadGraphS/aMCONLO package. The ti process was generated at tree-level, merged to 
one jet via the MLM method, and normalized to the NLO cross section. We show lu-equivalent errors, derived according to 
the Poisson distribution. For the case of bbZ in the ‘ideal’ parametrization, the ’<’ indicates the Icr-equivalent region, since no 
events were obtained after analysis. 


due to large backgrounds originating from the processes 
bbZj and 

The last final state considered was that containing 
(bb) and 2 leptons. This receives contributions from 
{bb){W~^W~), (56)(r+T“) and The former 

two channels possess large missing energy due to neutri¬ 
nos. The latter channel is not expected to be associated 
with significant missing energy, with the muons coming 
directly from the decay of the Higgs boson. The anal¬ 


ysis was split into two signal regions, depending on the 
amount of missing energy. The signal region targeting 
{bb){W^W~) and (66)(r+T“) with leptonic hnal states 
provides promising results, with few tens of events versus 
a few hundred background events at 3 ab“^ of integrated 
luminosity. This can be further improved by designing 
a method that captures the individual features of the 
two sub-channels, as well as using flavour information 
for the leptons. The other signal region, designed for the 





















{bb){fi~^ fi~) channel yields few events at 3 ab“^, with sig- 
nihcantly larger backgrounds. Potential enhancement of 
the latter channel could arise from an improvement in 
the resolution of the muon momenta that would allow a 
tighter mass window around their invariant mass. 

A boost in the importance of all the channels contain¬ 
ing a h —>■ 66 would similarly be obtained by improvement 
of the resolution of the 6 -jet momenta, as the mass win¬ 
dow around the Higgs boson mass could then be shrank 
down to beyond the 30-50 GeV region that has been con¬ 
sidered here. Additionally, in a final analysis, performed 
with a more complete FCC-hh detector design in mind, 
boosted decision tree or neural network methods would 
increase significances by taking into account the intri¬ 
cate correlations between the observables, going beyond 
the simple “rectangular” cuts that we have applied here. 
Moreover, a complete study should take into account re¬ 
alistic estimates of the theoretical uncertainties due to 
higher order corrections, as well as improved Monte Carlo 
descriptions for the generation of signal samples that will 
become available in the future. 

Finally, the initial investigation of the processes consid¬ 
ered here, along with the hh —)■ ( 66 ) ( 77 ) process, indicate 


that the study of Higgs boson pair production at a fu¬ 
ture hadron collider with 100 TeV centre-of-mass energy 
would greatly benefit with the collection of 10 or 30 ab“^ 
of integrated luminosity. 
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